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ABSTRACT 

We use images of high spatial, spectral and temporal resolution, obtained using both ground- and space- 
based instrumentation, to investigate the coupling between wave phenomena observed at numerous heights in 
the solar atmosphere. Analysis of 4170 A continuum images reveals small-scale umbral intensity enhance- 
ments, with diameters ^0'.'6, lasting in excess of 30 minutes. Intensity oscillations of «3 minutes are observed 
to encompass these photospheric structures, with power at least three orders-of-magnitude higher than the sur- 
rounding umbra. Simultaneous chromospheric velocity and intensity time series reveal an 87 ± 8° out-of-phase 
behavior, implying the presence of standing modes created as a result of partial wave reflection at the transition 
region boundary. We find a maximum wave guide inclination angle of «40° between photospheric and chromo- 
spheric heights, combined with a radial expansion factor of <76%. An average blue-shifted Doppler velocity 
of wl.5 kms -1 , in addition to a time lag between photospheric and chromospheric oscillatory phenomena, 
confirms the presence of upwardly-propagating slow-mode waves in the lower solar atmosphere. Propagat- 
ing oscillations in EUV intensity are detected in simultaneous coronal fan structures, with a periodicity of 
1 72 ± 1 7 s and a propagation velocity of 45 ± 7 km s "~ 1 . Numerical simulations reveal that the damping of the 
magneto-acoustic wave trains is dominated by thermal conduction. The coronal fans are seen to anchor into the 
photosphere in locations where large-amplitude umbral dot oscillations manifest. Derived kinetic temperature 
and emission measure time-series display prominent out-of-phase characteristics, and when combined with the 
previously established sub-sonic wave speeds, we conclude that the observed EUV waves are the coronal coun- 
terparts of the upwardly-propagating magneto-acoustic slow-modes detected in the lower solar atmosphere. 
Thus, for the first time, we reveal how the propagation of 3 minute magneto-acoustic waves in solar coronal 
structures is a direct result of amplitude enhancements occurring in photospheric umbral dots. 
Subject headings: MHD — Sun: chromosphere — Sun: corona — Sun: oscillations — Sun: photosphere — 
sunspots 



1. INTRODUCTION 

Early white-light eclipse photographs of the Sun revealed 
elongated, faint col umns of enhanced density stretc hing far 
out into the corona (Ivan de Hulsj[T950l fSaitol lT965h . These 
structures, now commonly referred to as coronal plumes, 
can be viewed over a wide range of wavelengths, i n par- 
ticular the extreme ultra violet (EUV; iBohlin et alJ 119751: 
I Ahmad & Withbroell 19771) . Plumes are just one type of EUV 
feature that are seen to exist in the solar corona. Other 
examples include coronal loop and fan structures, which 
are observed to outline the coronal magnetic-field topol- 
ogy, and demonstrate a wide range of oscillatory behaviour 
(e. g.. lAschwanden et al.ll 1 9991 120021: iNakariakov et al 1 1 1 999 1 : 
Jess et al.1 l2008ct lOfman & Wand l2008t IVerth et al.1 12008L 
Van Doorsselaere et al.ll2008HBallai et al.ll201 ll) . 

One of the first studies which uncovered propagat- 
ing wave phenome n a in coronal structures was that by 
iDeforest & Gurmanl (119981) . These authors utilised the 
Extreme-ultraviolet Imaging Telescope on board the Solar 
and Heliospheric Observatory spacecraft to identify quasi- 
periodic perturbations i n the brightness of 171 A images. 
iDe Moortel et al.l (1200 0) undertook a similar study using the 
higher spatial re solution Transition Region and Coronal Ex- 
plorer (TRACE: lHandv et al.lll999l) . and concluded that these 
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oscillations were signatures of slow magneto-acoustic waves, 
which propagate upwards along the coronal waveguides with 
velocities of 70-165 kms -1 and periods in the range 180- 
420 s. Energy estimates for these motions exhibit an incred- 
ibly wide range of values, typically 10 2 -10 5 erg s cm~ 2 s _1 
(IDeforest & Gurmanlll998l:lDe Moortel et al.ll2000l) . 

Since magnetic fields play an essential role in 
plume/fan/coronal-loop formation and structuring, they 
are often modelled using magnetohydrodynamic (MHD) 
equations (iDel Zanna et al.l 119971) . Utilising no n- linear , 
two-dimensional MHD simulations, Of man et al.l (119991) 
were able to replicate previous observational results, and 
concluded that outward-propagating slow magneto-acoustic 
waves may be able to contribute significantly to the heat- 
ing of the lower corona through compressive dissipation. 
Furthermore, theoretical modelling has suggested that the 
propagation characteristics of a magneto-acoustic intensity 
perturbation depends on a number of factors, including 
the dissipation of the wave ene rgy (lKlimchuketal.ll2004t 
IDe Moortel & N akariakovl 120121) . Numerical simulations 
indicate that thermal conduction may be the dominant damp- 
ing mechanism behind the dissipati on of magneto-a c oustic 
wave energy in the solar corona dOfman & Wand 120021: 
I De Moortel & Hoodl l200l I200H iMendoza-Bricefio et all 
12004ft . 

Through examination of coronal structures in the 
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FIG. 1 . — Simultaneous images of the lower solar atmosphere, obtained through 4170 A continuum (left), Ha core (middle), and Ca II 8542 A core (right) filters 
at 13:32 UT on 2011 July 13. The Ca II 8542 A image displays Doppler-compensated intensities, while the dashed white box outlines the ROSA/HARDcam 
field-of-view. Axes are in heliocentric arcseconds, where 1" ss 725 km on the solar surface. 



close proximity of a ctive regions, Fludra (2001) and 
lMarsh& WalshT (12006) were able to reveal how structures 
situated above sunspot regions displayed intensity oscil- 
lations with a period of the order of three minutes, while 
oscillations in "non- sunspot" structures dem onstrated much 
longer periodicities dDe Moortel et al.l 120021) . The authors 
concluded that the most likely explanation for the observed 
longitudinal waves revolves around a driver directly exciting 
the magnetic footpoints. This scenario requires the magneto- 
acoustic wave trains to be able to propagate from the lower 
solar atmosphere, through the transition region, and into the 
corona. Util ising numerical simulati o ns of the Sun's lower 
atmosphere, Kho menko & Co llados (2006) have revealed 
how longitudinal oscillations, generated with a periodicity 
of ^3 minutes in sunspot umbrae, can readily propagate 
upwards from the photosphere and into the chromosphere. 

Previously, three minute umbral oscillations have been 
notoriously di fficult to detect at photospheric heights. 
Balth asar et al.l (11987b were unable to detect any photospheric 
signatures using the Locarno solar s t ation at the Gottingen 
Observatory, while Lites & Thomas (1985) suggested that 
they may get lo st in the noise, as a result of their very low am- 
plitude. Indeed, Nagashi ma et al.l d2007 ) utilised Hinode/SOT 
image sequences to show how oscillatory power, at all fre- 
quencies , is significantly reduced in sunspot umbrae. More 
recently, Koban ov et al.l (120081 1201 11) have not only detected 
photospheric three minute oscillations, but the authors also 
claim that the location of maximum chromospheric power 
also corresponds to a co-spatial decrease in power of the pho- 
tospheric osdllatiojis.How^er ! _Jhe spatial resolution ob- 
tained by Koba nov et all ( 120081 l20Tlh was on the order of 1", 
so precise diagnostics of the exact umbral structures display- 
ing three minute periodicities was impossible. 

In this paper, we utilise ground- and space-based instru- 
mentation, with high spatial, temporal and spectral resolu- 
tion, to investigate the origin of 3 minute magneto-acoustic 
waves observed in EUV images of coronal fan structures. 
We employ a multi-wavelength approach to study the photo- 
spheric counterpart of these coronal phenomena, and analyse 
the resulting wave propagation characteristics from the pho- 
tosphere, through the chromosphere, and out into the corona. 



2. OBSERVATIONS & DATA PREPARATION 
2.1. Ground-based Data 

The observational data presented here are part of a se- 
quence obtained during 13:32 - 14:04 UT on 2011 July 13, 
with the Dunn Solar Telescope (DST) at Sacramento Peak, 
New Mexico. We employe d the Rap i d Os cillations in the 
Solar Atmosphere (ROSA; iJess et all 1201 Obi) camera sys- 
tem to image a 69" x 69" region encompassing active re- 
gion NOAA 11250, positioned at heliocentric co-ordinates 
(-146", -486"), or S26.3E10.1 in the conventional north- 
south-east-west co-ordinate system. A spatial sampling of 
'.'069 per pixel was used for the ROSA cameras, to match 
the telescope's diffraction-limited resolution in the blue con- 
tinuum to that of the CCD. This results in images obtained 
at longer wavelengths being slightly oversampled. However, 
this was deemed desirable to keep the dimensions of the field- 
of-view the same for all ROSA cameras. 

A recent addition to the DST's imaging capabilities is a 
new high quantum efficiency device, the Hydrogen-Alpha 
Rapid Dynamics camera (HARDcam). This camera, an iXon 
X3 DU-887-BVQ model manufactured by Andor Technol- 
ogy, consists of a back-illuminated, 512x512 pixel 2 electron- 
multiplying CCD, with a quantum efficiency exceeding 95% 
at 6500 A. As a result, this camera is best suited to imag- 
ing in the red portion of the optical spectrum. The triggering 
and readout architectures are identical to the first-generation 
ROSA instrumentation, allowing HARDcam to be seamlessly 
integrated with the existing setup. To optimise the exploita- 
tion of its high quantum efficiency, HARDcam was employed 
behind a 0.25 A Ha core filter, incorporating a spatial sam- 
pling of 0'.'138 per pixel, providing a field-of-view size (71" 
x 71") comparable to existing ROSA image sequences. 

In addition to ROSA and HARDcam observatio ns, the In- 
terfer ometric BIdimensional Spectrometer (IBIS; ICavallinil 
2006) was used to simultaneously sample the Ca II absorp- 
tion profile at 8542.12 A. IBIS employed a spatial sampling 
of 0'.'097 per pixel, which allowed ROSA's near square field- 
of-view to be contained within the circular aperture provided 
by IBIS. Thirteen discreet wavelength steps, with ten expo- 
sures per step to assist image reconstruction, were used to 
provide a complete scan cadence of 43.4 s. A whitelight cam- 
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FIG. 2. — Top Left: An AIA 171 A image, taken at 13:32 UT, including a dashed white box to indicate the ROSA/HARDcam field-of-view. Top Right: An 
AIA 4500 A context image, interlaced with a ROSA 4170 A continuum snapshot, following sub-pixel co-alignment. Bottom Panels: Excellent co-alignment 
between lower and upper solar atmosphere datasets allow for precise spatial diagnostics, as demonstrated by the interlacing of 4170 A continuum (left) and Ha 
core (right) images, with the coronal emission seen in the AIA 171 A bandpass. 



TABLE 1 

ROSA/HARDC AM/IBIS FILTER AND CADENCE OVERVIEW. 

Filter Used Exposure Time Frames per Reconstructed 

(ms) Second Cadence (s) 

Continuum (4170 A) 7 303 1.06 

Ha core 35 27.9 1.26 

Call (8542.12 A) 250 a 3.0 6 43.4 

a : Average frames per second including Fabry-Perot tuning time 
b : Cadence of a complete Ca II profile scan 

era, synchronised with the IBIS feed, was utilised to assist the 
processing of narrowband images. Full details of the obser- 
vations presented here, including filters and exposure times 
used, can be found in Table Q] while sample images can be 
viewed in Figure Q] 

During the observations, high-order adaptive optics 
(Rimmele 2004) were used to correct wavefront defor- 
mations in real-time. The acquired images were fur- 
ther improved thr ough speckle reconstruction algorithms 
dWoger et alj|2008l) . utilizing 32 ->■ 1 restorations for the G- 
band and 4170 A continuum images. The remaining HARD- 
cam and IBIS images were processed with 35 —> 1 and 
10 ^ 1 restorations, respectively. Post-reconstruction ca- 
dences are displayed in the fourth column of Table Q] A full 



image-reconstructed IBIS scan through the Ca II absorption 
line includes a blueshift co rrection, required d ue to the use of 
classical etalon mountings dCauzzi et al.ll2.Q08h . Atmospheric 
seeing conditions remained excellent throughout the time se- 
ries. However, to insure accurate co-alignment in all band- 
passes, broadband time series were Fourier co-registered and 
de-stretched using a 40 x 40 g rid, equating to a «1_^7 separa- 
tion between spatial samples dJess et al.l 2007a, 2008a). Nar- 
rowband images, including those from IBIS, were corrected 
by applying destretching vectors established from simultane- 
ous broadband reference imag es (Reard on & Cavalfin"ill2008l : 
Uess et alJl2010llReardon et al.ll201 lh . " 

2.2. Space-based Data 

The Atmospheric Imaging Assembly (AIA; iLemen et alJ 
1 201 lh onboard th e Solar Dynamics Observatory (SDO; 
Pesne llet alj201 II ) was utilised to provide simultaneous EUV 
images of active region NOAA 11250. This instrument im- 
ages the entire solar disk in 10 different channels, incorpo- 
rating a two-pixel spatial resolution of 1 '.'2, and a cadence of 
12 s. Here, we selected 7 EUV datasets spanning 13:30 - 
14:05 UT on 2011 July 13, consisting of 175 images in each 
of the 94A, 131 A, 171 A, 193 A, 211 A, 304A, and 335A 
channels. In addition, one contextual 4500 A continuum im- 
age, acquired at 14:00 UT, was obtained for the purposes of 
co-aligning AIA datasets with images of the lower solar at- 
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FIG. 3. — Left: A time-averaged 4170 A continuum image, where a solid green contour defines the outline of the sunspot umbra. Middle and Right: A time- 
averaged close-up of the sunspot umbra, following isolation from the surrounding penumbral and granulation structures. Locations of the umbral dots analysed 
here are indicated by red "UDx" markings, where the "x" refers to a number. Green, red, and blue contours indicate the locations where photospheric intensities, 
chromospheric Hct/Ca II line core intensities, and Doppler velocities, respectively, display oscillatory power exceeding 1000 times the background umbral value, 
for a periodicity of approximately 3 minutes. 



mosphere. 

The EUV bandpasses were specifically chosen to cover 
a wide range of transition region and coronal temperatures, 
spanning approximately 50 000 K - 7 MK, under non-flaring 
conditions. Transition region imaging is covered by the He II- 
dominated 304 A bandpass, with a typical formation tem- 
perature of ^50 000 K. The selected coronal channels, con- 
sisting of the 94A, 131 A, 171 A, 193 A, 211 A, and 335A 
bandpasses, demonstrate typical effective temperatures of ap- 
proximately 7.0 MK, 0.4 MK, 0.7 MK 1-6 MK, 2.0 MK 
and 2 .8 MK respectively dO'Dwver et al.ll2010UBrooks et al] 
1201 ll) . Thus, the AlA datasets, in conjunction with our 
ground-based observations of the lower solar atmosphere, 
provide us with the ideal opportunity to investigate the cou- 
pling of lower atmospheric (photosphere and chromosphere) 
phenomena with their multi-million degree coronal counter- 
parts at unprecedented spatial and temporal resolutions. 

The AIA data were processed using the standard 
aia_prep routine, and include the removal of energetic par- 
ticle hits, in addition to the co-registration of images from 
different wavelengths on to a common plate scale. Sub- 
sequently, 200" x 200" sub-fields were extracted from the 
processed data, with a central pointing close to that of the 
ground-based image sequences. A sub-field image, includ- 
ing an outline of the field-of-view obtained with ground- 
based observations, is shown in the upper-left panel of Fig- 
ure |2] Using the AIA 4500 A context image to define abso- 
lute solar co-ordinates, our ground-based observations were 
subjected to cross-correlation techniques to provide sub-pixel 
co-alignment accuracy between the imaging sequences. To 
do this, the plate scales of our ground-based observations 
were first degraded to match that of the AIA imag efl Next, 
squared mean absolute deviations were calculated between 
the datasets, with the ground-based images subsequently 
shifted to best align with the AIA reference image (see the 
lower panels of Figure|2]for interlaced examples). Following 
co-alignment, the maximum x- and y-displacements are both 
less than one tenth of an AIA pixel, or 0'.'06 («45 km). 

3. DATA ANALYSIS AND DISCUSSION 

2 Data analysis was performed on full-resolution (i.e. non-degraded) im- 
age sequences 



Using a combination of ground- and space-based datasets, 
we were able to analyse image sequences from the photo- 
sphere through to the corona. Due to the wide range of at- 
mospheric heights covered, this section is divided into sub- 
sections in which we discuss different regions of the solar at- 
mosphere. 

3.1. Photosphere 

To investigate oscillatory phenomena occurring in the 
sunspot umbra, we first masked out all other areas of the 
field-of-view. This was deemed desirable so that intrinsically 
faint umbral structures would become much more apparent in 
scaled intensity images. To isolate the umbra, we first had 
to create a binary map, whereby umbral pixels were given a 
value of '1', and non-umbral pixels assigned a value of '0'. 
This form of binary map was created by first averaging the 
4170 A continuum images over the entire 32 minute duration 
of the dataset. Next, the umbral pixels were defined as those 
with an intensity below 45% of the median granulation in- 
tensity, and subsequently assigned a value of '1'. This ac- 
curately defined the perimeter of the dark umbra. However, 
some small-scale structures, which existed inside the umbra 
(e.g., umbral dots; UDs), exhibited higher relative intensity 
values. Therefore, to include all structures which existed in- 
side the umbral boundary, all pixels which lay inside the outer 
perimeter were assigned a value of '1', while all other pixels 
were given a value of '0'. The time-averaged 4170 A contin- 
uum image, including the outline of the umbral binary map, 
is shown in the left panel of Figure [3] 

Examination of the time-averaged umbral intensity map 
(middle panel of Figure [3} revealed a collection of near- 
circular brightenings. These brightenings are consistent with 
the signatures of UDs, whereby they typically exhibit diam- 
eters between 0'/2 - 0'.'8, and intensities 1 3 - 1.5 times 
brighter than the backgrou n d umbral value dDenkei 119981 ; 
iTritschler & Schmidtl 12001: iSobotka & Hanslmeieil l2Q05h . 
Due to their sharp appearance in a time-averaged intensity 
image, it hints at the persistent nature of such structures in 
the same spatial loc ation, either through l ong lifetimes (e.g., 
30 minutes or more; Sobotka et al. 1997113), or by the reap- 
pearance of new UDs in the same spatial location. Analysing 
simultaneous G-band and continuum filtergrams, Rimmele 
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FIG. 4. — Upper: Intensity time series of the four umbral dots labelled UD1 
(red line), UD2 (blue line), UD3 (green line), and UD4 (pink line) in FigurefJ] 
following temporal filtering using a 2 - 4 minute bandpass filter. The black 
line displays the intensity of the umbral region away from UD structures, 
devoid of long-lived or transient brightenings, and scaled (by a factor of 5) 
for clarity. Lower: The relative phase angles between the 3 minute umbral 
dot and background umbral oscillations, plotted as a function of time. The 
red, blue, green, and pink lines represent umbral dots UD1, UD2, UD3, and 
UD4, respectively, while a black dashed line highlights a phase angle of 0° . 

(2008) found that UDs, when observed in the G-band, dis- 
played dark fine-structuring towards their core. Contrarily, 
this author found that these dark features are not as readily 
visible in the continuum, and concluded this may be due to 
the reduced resolution of the longer wavelength and longer 
exposure images. However, our continuum images were ob- 
tained at a shorter wavelength than the G-band (4170 A versus 
4305 A), in addition to being captured with a much shorter ex- 
posure time (7 ms versus 20 ms). Thus, we cannot attribute 
a lack of dark fine-scale structuring in the UDs observed in 
the 4170 A continuum to either reduced spatial resolution, 
or exp osure-time related smearing of the images. iJess et all 
(120121) recently determined the contribution functions of the 
G-band and 4170 A continuum filter bandpases used at the 
DST, and concluded that the blue continuum is formed at a 
height of ^25 km, while the G-band is formed approximately 
100 km higher. This atmospheric height separation between 
the two bandpasses, in addition to the appearance of diatomic 
CH molecules in the G-band images, may contribute to the 
presence of dark cores observed in G-band UDs. 

Since we are concerned with the larger macro-scale fluc- 
tuations in umbral intensity, rather than the very small-scale 
structuring (e.g., dark features within UDs) which may be 
close to, if not overlap with the telescope's diffraction limit, 
our photospheric umbral time series was created purely from 
consecutive 4170 A continuum images. These images were 
multiplied by the umbral binary map, producing a dataset 
which could be anal ysed with Fourier techniques. Following 
the methodolog y of | Jess et al„ (2007a b), the wavelet analy- 
sis routines of lTorrence & Compol (119981) were applied to the 
4170 A continuum umbral time series to search for the pres- 
ence of oscillatory behaviour. Considerable oscillatory power 



was present throughout the sunspot umbra, with dominant 
peaks in the Fourier spectrum at 3 and 5 minutes. This is con- 
sistent with the generalised p-mode spectrum of sunspot oscil- 
lation s dLites etal.ll 19821: lLitesl[l984 Il986t iBogdan & Judgel 
2006). However, when locations of high oscillatory power 
were examined, it became evident that the UDs displayed sig- 
nificant enhancements of wave amplitude, at periodicities of 
approximately 3 minutes. The middle panel of Figure[3]shows 
the spatial locations where oscillatory power at approximately 

3 minutes is in excess of 1000 times the background um- 
bral value. To test whether these long-lived umbral bright- 
enings oscillate in phase with the surrounding umbra, tempo- 
rally filtered time series of these structures were first created 
by passing each spatially-averaged UD intensity through a 2 - 

4 minute bandpass filter. The resulting lightcurves, dominated 
by power corresponding to the 3 minute p-mode oscillations, 
are displayed in the upper panel of Figure |U along with a fil- 
tered time series representing a region of the quiet umbra de- 
void of any transient brightenings. It is clear that not only do 
the peaks and troughs of the UD oscillations appear similar in 
time, but these signatures also closely follow the oscillations 
originating from within the background umbra. 

To quantify any small-scale differences in the oscillating 
UD time series, a phase difference analysis with respect to 
the background umbra was performed. The lower panel of 
Figure|4]reveals the phase difference between each of the four 
UD time series, and an isolated region of the background um- 
bra (i.e. the solid black line in the upper panel of Figure [H. 
It is clear that a preference exists for the UDs to oscillate in 
phase with the background umbra. Maximum deviations in 
the phase angle reach approximately ±40°, although these 
may be due to slight drifts in oscillation period between the 
respective time series. For example, the phase angle related to 
UD 1 (red line in the lower panel of Figure HI drifts a total of 
«80° over a 700 s time duration, providing an average shift of 
«0.11° s _1 . This drift can be explained by a period discrep- 
ancy of only 20 s between the UD and the background umbra, 
something which is well w ithin the normal range of p-mode 
frequencies (Thomas 1985). Thus, we interpret the umbra as 
a single oscillating "drum skin", which not only induces wave 
motion in its inherent UDs, but also causes these structures to 
oscillate in phase wi t h one another. 

While iChou et al.1 (120091) and iStangalini et all J20TI have 
shown that 3 minute magneto-acoustic power is suppressed in 
sunspot umbrae as a result of local absorption and emissiv- 
ity reduction, their spatial resolution was insufficient to allow 
studies of the smallest sunspo t features (e.g., umbral dots). 
Recently, Shelyag et al. (2009) have shown that on small 
spatial scales, the curvature and strength of magnetic field 
lines have substantial influence on the efficiency of magneto- 
acoustic wave propagation. Indeed, these authors suggest 
that wave power can actually be amplified under certain at- 
mospheric conditions. Following on from this, our observa- 
tional data suggests that UDs are able to enhance the back- 
ground p-mode power by at least three orders-of-magnitude 
(see e.g., middle panel of Figure|3}. This may be a result of the 
magnetic field lines acting as effic ient conduits for magneto- 
acoustic wave propagation (e.g.. ISinghlll992t IShelvag et alJ 
l2006tlKhomenko et al. 2008; Erdelyi & Fedun 2013). 
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3.2. Chromosphere 

Chromospheric information comes from the Ha core imag- 
ing dataset, in addition to the Ca II 8542 A spectral imag- 
ing scans taken with IBIS. Furthermore, Doppler shifts of the 
Ca II profile minimum allo w a series of two- dimensional ve- 
locity maps to be generated dJess et alj20 10a). The Ca II core 
image shown in the right panel of Figure Q] is a true intensity 
map, created by establishing the line-profile minimum at each 
pixel. By displaying Doppler-compensated line-centre inten- 
sities, rather than rest-wavelength intensities, brightness vari- 
ations throughout the image are more indicative of the source 
functi on than of the veloc ities present in the line-forming re- 
gion (Leenaarts et al. 20lO). 

In a process identical to that applied to the photospheric 
image sequence, the Ha and Doppler-compensated Ca II data 
sets were multiplied by the umbral binary map, and subse- 
quently analysed using wavelet techniques. As for the 4170 A 
continuum time series, considerable oscillatory power was 
present throughout the entire sunspot umbra, with dominant 
peaks in the Fourier spectru m at 3 a n d 5 m inutes. However, 
remaining consistent with iThomasi ([1985), chromospheric 
power at periodicities of ~5 minutes was much reduced when 
compared to those found in the photospheric umbra. To inves- 
tigate whether oscillatory power detected in the chromosphere 
can be related to similar phenomena found in the photosphere, 
the locations of high oscillatory power were examined. The 
right panel of Figure [3] outlines the spatial locations where 
power, in both intensity and velocity signals at approximately 
3 minutes, is in excess of 1000 times the background umbral 
value. Similar to the locations of highphotospheric power dis- 
played in the middle panel of Figure|3] several distinct groups 
manifest in a horseshoe shape in the south-west quadrant of 
the sunspot umbra. 

The locations of significant chromospheric power have 
two distinct differences when compared to their photospheric 
counterparts. First, their spatial positions are slightly offset 
from those displayed in the middle panel of Figure [3] Sec- 
ondly, the spatial sizes of high chromospheric power are sub- 
stantially larger than those found in the photosphere. These 
effects can be explained by the geometry of the magnetic field 
lines which stretch outwards from the solar surface. For ex- 
ample, high photospheric power encompassing UD2 (mid- 
dle panel of Figure [3} has a local maximum at the helio- 
centric co-ordinate (— 133'.'7, — 483'.'7), while the chromo- 
spheric local maximum is at (— 133 '.'1, —482 '.'2). Thus, the 
offset between these two maxima is wl'/6, or «1100 km. A 
separation between the 4 170 A continu um and the Ha/Ca II 
core formation heights of ^1800 km (IVernazza et al.1 Il98lt 
Uess etalf2 012) requires an inclination angle of w40° to sup- 
port the assumption that the oscillatory sign als are from a 
contin uation of the same solar structure (e.g., ICenteno et al.l 
2006). This is the largest offset present in our observations, 
and fo rms our upper limit of the magnetic field inclination 
angle. IRempel et al.l (120091) have recently shown that the in- 
clination angles (to the vertical) of magnetic flux tubes within 
sunspot umbrae can be quite large (>45°), with these struc- 
tures often be coming horizontal n ear the penumbral boundary. 
Furthermore, Mars h et al.l d2009l) utilized stereoscopic obser- 
vations of propagating slow-mode waves in coronal structures 
to infer an absolute inclination angle ^40°. Thus, an inclina- 
tion angle of <40° supports the interpretation that oscillatory 
behaviour detected at photospheric and chromospheric layers 
are directly related by the magnetic field lines which extend 
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FIG. 5. — Upper: An unfiltered time series, created by spatially averag- 
ing Ca II intensity (solid line) and velocity (dashed line) signals, along with 
Ho core intensities (red line), originating from within the contours outlining 
'UD2' in the right panel of FigurefJ Middle: The same time series from the 
upper panel, following temporal filtering through a 2 - 4 minute bandpass fil- 
ter. A —90° phase delay between velocity and intensity signals indicates the 
presence of oscillations which are magneto-acoustic in nature, while a close 
agreement between Ha and Ca II intensities is a result of their similar for- 
mation heights. Lower: The relative phase angles between the photospheric 
and chromospheric umbral dot oscillations, plotted as a function of time. The 
black, red, green, and blue lines represent umbral dots UD1, UD2, UD3, and 
UD4, respectively, while a black dashed line highlights a phase angle of 0° . 

upwards from the solar surface. 

An increase in the area of the oscillating regions can be 
associated with a physical expansion of the magnetic flux 
tubes as a function of atmospheric height. Through exam- 
ination^ intens e magnetic field concentrations (>1000 G), 
Jess et al. (2009) were able to show that the radial expansion 
of magnetic flux tubes between photospheric and chromo- 
spheric heights can be as large as a factor of 2. Assuming 
these magnetic flux tubes demonstrate an expanding cylin- 
drical geometry, doubling the radial dimension will result 
in an area increase of 400%. Continuing with the previous 
example, the region of high oscillatory power encompass- 
ing UD2 covers 326 pixels (0.8 Mm 2 at 2500 km 2 pixel -1 ) 
and 256 pixels (2.5 Mm 2 at 10 000 km 2 pixel -1 ) at pho- 
tospheric and chromospheric heights, respectively. As a 
result, an expansion of only 310% is observed, meaning 
our observations are well within the limits of previous ex- 
panding magnet ic flux tube models (e.g., iMein et aD 120071: 
Ruderman et alj|2008t IShelyag et al.l l2mot iKarami & Bahari 
201lHFedunetal.ll2011l) . " 

The specific mode of oscillation can be determined through 
investigation of the coupling between intensity and ve- 
locit y signals. Here, we adopt the V — I conven- 
tion dDeubner & Flecklll989l:lFleck & Deubnerll 19891) . which 
shows the delay of maximum intensity (/) with respect to 
maximum blueshift velocity (V). Thus, a wave which has 
a velocity signal trailing its intensity signal by 1/4 of a pe- 
riod will have a V — I phase angle of —90°. Due to the 
Ca II absorption line being sensitive to temperature fluctua- 
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tions (Bee be&Johnson| [T969). we follow the common prac- 
tice to a dopt the inte nsity (/) as a proxy for the local tem- 
perature. iMeinl d 19771) has shown that when waves propagate 
along moderately inclined flux tubes, the phase lag between 
fluctuations in velocity and temperature are the same as those 
found in purely acoustic waves which have been modified by 
gravity. Therefore, in a fully adiabatic scenario, theory pre- 
dicts a V — I phase angle approaching 0° for running acoustic 
waves, which can increase to —90° whe n aspects of wave re- 
flection create standing acoustic modes (|Hofmann et alJ 19961 : 
IA1 et all 1 9981 iNi gam & Kosovichevll999l) . Under isothermal 
conditions, the phase lag can fu rther increase up to —180° 
(Miha las & Toomrd[l98ll 119821) . Examination of the upper 
panel in Figure [5] reveals how the unfiltered velocity time 
series (dashed line) trails the co-spatial intensity lightcurve 
(solid line) by approximately 1/4 of a wave period. This ef- 
fect becomes even more pronounced when the same time se- 
ries are passed through a 2 - 4 minute bandpass filter (mid- 
dle panel of Figure |5). Chromospheric waves detected within 
the contours of the right panel of Figure [3] have a spatially 
and temporally averaged phase angle of —87 ± 8°, suggest- 
ing these waves are magneto-acoustic in nature, with char- 
acteristics consis tent with standing acoustic modes (iDeubnerl 
ll974tlCramlll978l) . The generation of a chromospheric stand- 
ing wave may be the result of a portion of the acoustic wave 
energy being reflected back at the transition region boundary 
( Schm itz & Fleck1[l99l iNakariakov etld1l2~004t iFedun et al.1 
120111) . 

The spatial and temporal averaging of velocity signals 
found in the locations of high chromospheric power results 
in a net blueshift velocity of wl.5 kms -1 (upper panel of 
Figure |5). While aV — I phase angle of —87 ± 8° implies 
the presence of wave reflection, a net blueshift velocity may 
also indicate that a portion of the magneto-acoustic wave is 
propagating in the upward direction. To examine the prop- 
agation characteristics of the magneto-acoustic waves asso- 
ciated with the UDs, spatially-averaged Ha intensity time- 
series were constructed for each of the regions contoured in 
the right panel of Figure [3] and subsequently passed through 
a 2 - 4 minute bandpass filter to isolate oscillatory phenomena 
around 3 minutes. A phase difference analysis was performed 
between these chromospheric Ha lightcurves and their cor- 
responding pho tosp heric 4170 A continuum counterparts cre- 
ated in Section 13.11 and displayed in the upper panel of Fig- 
ure |4] A preference for negative phase angles is shown in 
the lower panel of Figure [5] which indicates the oscillatory 
signals are first observed in the photospheric 4170 A contin- 
uum, before later being detected in the chromospheric Ha 
time-series. By averaging the derived phase angles over all 
UDs, an average phase lag of —43 ± 12° is found, imply- 
ing the presence of upwardly propagating waves. Using the 
dominant periodicity of «3 minutes, a relative phase angle of 
—43° equates to a time lag of w22 s. This is consistent with 
the photos phere-to-chromosphe re time-lag measurements de- 
tailed in Kobanovetal. (2011). However, due to a height 
separation of ^1800 km betw een the 4170 A cont i nuum and 
Ha c ore image sequences (IVernazza et al.1 11981b iJess et all 
12012b . a time lag of «22 s requires a phase speed exceed- 
ing 80 kms -1 . This is clearly unfeasible since a velocity 
this large will be supersonic, and violates our interpretation 
that the observed waves are magneto-acoustic slow modes. 
However, because the wave trains exist at the start of our ob- 
serving sequence, and continue to be observed as the time- 



series finishes, a factor of n2ir (n360°), where n is an integer, 
may be absent from our derived phase angle. By considering 
n = 1, the phase lag increases to «— 403°, with the result- 
ing phase speed reducing considerably to «8 kms -1 . Only 
by observing the start and/or end of a propagating wave train 
can the exact time lag (and associated phase speed) be conclu- 
sively determined. Nevertheless, the detected waves are best 
described as upwardly-propagating magneto-acoustic modes, 
which travel along magnetic flux tubes inclined to the verti- 
cal by less than ft!40 o . Such an i nclination angle may explain 
whv lKobanov et all (l2008l 120111) were unable to correlate (on 
a pixel-by-pixel basis) high chromospheric oscillatory power 
to that occurring in the underlying photosphere. 

Our Ha dataset has a much higher cadence than the Ca II 
core image sequence (1.26 s instead of 43.4 s). Thus, exami- 
nation of the Ha intensity time series can substantially reduce 
the associated errors when deriving the periodicities of the 
propagating magneto-acoustic waves. As can be seen in Fig- 
ure [5] the Ha data displays prominent intensity oscillations, 
most of which last for the entire duration of the time series. A 
spatial and temporal averaging over the regions encompassed 
by the red contours in the right panel of Figure [3] yields a 
peak chromospheric periodicity of 168 ± 7 s, consistent with 
the generali sed jp-mode spectrum of chromospheric sunspo t 
oscillations dO'Shea et al.ll200Tl l2002t iBaneriee et al1l2002l) . 
The close resemblance between Ha and Ca II intensity time 
series suggests their formation heights are remarkably similar. 
This is also apparent by the similarities present in simultane- 
ous snapshots through their respective filters (e.g., FigureHJ. 

3.3. Transition Region and Corona 

Due to the reduced spatial resolution of AIA images (2 pixel 
resolution «1'.'2 «870 km), when compared with our simul- 
taneous datasets of the lower solar atmosphere (2 pixel reso- 
lution «0'/138 «100 km), it is imperative to concentrate on 
the larger-scale structures which will be apparent in all image 
sequences. Following inspection of the co-aligned images dis- 
played in the lower panels of Figure|2] it is clear that a number 
of fan structures rise out of the photospheric sunspot umbra, 
and into the corona. These structures are particularly visible 
as intensity enhancements in the coronal AIA images, where 
the 131 A, 171 A, 193 A, 211 A, and 335A bandpasses show 
similar structuring extending outwards in the south-west di- 
rection. However, these fan structures are mostly absent in 
the transition region 304 A images. 

To investigate whether an absence of the fan brightenings 
in the 304 A images is a consequence of the structures ly- 
ing outside of the filter's temperature response curve, differ- 
ential emission measure (DE M) techniques were employed 
(e.g.. lHannah & Kont ar 2012). Utilising the six coronal EUV 
channels on AIA, we were able to construct emission mea- 
sure (EM) and temperature (T e ) estimates of active region 
NOAA 11250, including its im mediate vicinity. Follow ing the 
methodology presented by lAschwanden et al] (1201 ll) . each 
batch of near-simultaneous EUV exposures allowed us to con- 
struct a time series of EM and T e variables, thus allowing 
the evolutionary changes in each sequence to be studied. Im- 
portantly, the temperature of the coronal fan structures cur- 
rently under investigation is in the range of 0.5 - 1.2 MK 
(right panel of Figure |6j. Their relatively cool temperature, 
when compared to the multi-million degree values found at 
the centre of the active region, probably manifest as a result 
of open magnetic field lines which cannot trap heated plasma. 



8 



D.B. Jess et al. 



Emission K/easure 



Temperature 



-400 



-450 



-500 



-550 




-400 



-200 -150 
Solar x (arcsecs) 




FIG. 6. — Two-dimensional emission measure (in units of cm -5 K -1 ; left) and temperature maps (in units of log(T e ); right), derived from near-simultaneous 
EUV snapshots of NOAA AR 11250. The coronal structures currently under investigation are located towards the south-west of the active region, where 
temperatures in the range of 0.6 - 1.2 MK dominate. 



An alternative explanation could revolve around their connec- 
tion with quasi-separatrix layers, which may be subject to a 
peculiar heating function (Schrijv er et al.ll2010l) . While the 
304 A bandpass has two distinct temperature response func- 
tions covering approximately 50 000 K and 1.5 MK, the fact 
that active region NOAA 11250 was positioned close to so- 
lar disk centre suggests that the resulting images should be 
dominated by the He II emission formed at ^50 000 K, with 
contributio ns from the 1 .5 MK Si XI 202.22 A emission line 
minimal (lO'Dwyer et aLllToiOh . This helps to explain why 
features near one million degrees (e.g., the fan structures ex- 
tending outwards from the underlying sunspot), are not read- 
ily apparent in the transition region data. 

Examination of time-lapse movies of coronal EUV images 
revealed clear and distinctly periodic outflows along the coro- 
nal fans. To quantify the associated periodicities and flow ve- 
locities, a series of one-dimensional slits were placed along 
the motion path in each of the coronal EUV channels. The re- 
sulting time-distance cuts reveal numerous propagating wave- 
fronts, indicated by straight, diagonal trends in the bottom 
panel of Figure [7j The AIA EUV images, in the temperature 
range of 0.4- 2.8 MK (incorporating the 131 A, 171 A, 193 A, 
21 1 A, and 335 A bandpasses), display a dominant periodicity 
of 172±17s, and a propagation velocity of 45±7 km s" 1 . Red 
and blue contours in the lower panel of Figure [7j outline the 
intensity signals present in the higher temperature 211 A and 
335 A bandpasses, respectively. It is clear that not only do the 
periodicities of the wave fronts closely resemble one another 
in both space and time, but the intensity gradients (and hence 
wave speeds) are also similar. Propagating wave fronts are 
observed co-spatially, and simultaneously in all but the 94 A 
AIA bandpass. The signal-to-noise of the 94 A channel is too 
low to extract a time series of sufficient quality for analysis. 
This may be compounded by the fact that the fan structures 
demonstrate a temperature much lower than the 94 A chan- 



nel's peak response. 

A periodicity of «172 s and a propaga t ion v elocity 
«45 kms -1 are consistent with iKing et all ([2003), who 
utilised the 171 A and 195 A filters onboard the TRACE satel- 
lite to reveal how the propagation velocity of such waves will 
not exceed the local sound speed, even with a large inclina- 
tion angle of the wave guides away from the observer's line- 
of-sight. The typical sound speed associated with the domi- 
nant Fe IX emission from the 171 A bandpass is ?»150 km s _1 
dDe Mo ortel 2006), which requires the fans currently under 
investigation to have an inclination angle exceeding 70° be- 
fore the observed wave motion would become supersonic. 
An inclination angle this large is highly unrealistic due to 
the location of the active region on the solar disk, in addi- 
tion to previous surveys of coronal magnetic geometries (e.g., 
lAschwanden et al.ll2002ll2008l 120091) . Thus, we can conclude 
that the observed wave phenomena is best described as prop- 
agating slow-mode waves. 

Since propagating wave fronts are detected in many of 
the AIA EUV bandpasses, emission measure and tempera- 
ture maps were subsequently investigated for the presence 
of similar signatures. As the emission measure maps cre- 
ated here are the sum of (squared) electron densities along 
a given line-of-sight, the total mass distribution, regardless of 
the local temperature, ca n be studied as a function of time 
( Aschwan den et all 1201 1). Furthermore, the derived pixel- 
by-pixel temperatures correspond to the peak of the DEM 
distribution for a given line-of-sight, thus allowing the en- 
tire tem perature range of the coronal plasma to be easily 
studied (lAschwanden & Boernerll2~01 ll) . Applying the same 
techniques used on the EUV image sequences, time-distance 
diagrams of the EM and T e variables derived for the fan 
structures reveal identical periodicities and propagation ve- 
locities as found in the EUV images (middle and upper pan- 
els of Figure |8). The emission measure is usually defined as 
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FIG. 7. — An AIA 171 A image (top), co-spatial with our field-of-view of 
the lower solar atmosphere. Dashed white lines outline a coronal fan struc- 
ture, originating from within the underlying sunspot umbra, where propa- 
gating wave phenomena are readily apparent. The perimeter of the sunspot 
umbra is highlighted by a solid red line, demonstrating how the coronal fan 
is anchored into regions of high chromospheric oscillatory power (see e.g., 
the right panel in Figure [5). A sample time-distance diagram of this fan is 
displayed in the lower panel, where km represents the sunspot umbra. Red 
and blue contours outline 211 A and 335 A intensities, respectively, which 
are 60% above the local quiescent background. A black dashed line repre- 
sents the propagation of a typical wave front, where the gradient provides the 
wave speed, typically 45 ± 7 km s — 1 . Wave phenomena which are co-spatial, 
co-temporal, and propagating with the same wave speed are readily apparent 
over a range of coronal temperatures. 

dKingston et alJ[l982l:lDovle et alj[l985l) . 

EM(T e ) = I nldh. 



(1) 



where n e is the electron density and h is an emitting depth 
along the observer's line-of-sight. Since the compressive 
phase of a magneto-acoustic wave mode will cause the emit- 
ting volume to decrease, thus reducing the overall emission 
measure, one would expect the EM time series to oscillate in 
phase with the corresponding EUV intensities. However, as 
the plasma becomes compressed, the associated kinetic tem- 
perature will increase, thus causing the T e signal to oscillate 
out-of -phase with the EM and EUV intensity signals. Indeed, 
perturbations in the EUV intensity, depicted in Figure [7j are 
found to occur in phase with the derived emission measure 
fluctuations. Furthermore, the lower panel in Figure[8]reveals 
how oscillations in T e are found to be 180° out-of-phase with 
respect to those detected in the EM signal. Quantitatively, the 
detected fluctuations in the EM and T e values are typically in 
the range 21.50±0.03 cm" 5 K" 1 and 0.55±0.03 MK, respec- 
tively. This strengthens our interpretation that the observed 
wave phenomena are best described as magneto-acoustic slow 
modes. 

3.4. Numerical Simulations 

To model our observational findings, we performed a 
series of ID numeri cal simulations us ing the Lagrangian- 
Remap code (LARE; lArber et all 1200 lb . in which thermal 
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FIG. 8. — Time-distance diagrams of the emission measure (in units of 
cm -5 K — 1 ; top) and temperature (in units of log(T e ); middle), derived for 
the fan structure where propagating wave phenomena in EUV intensity im- 
ages was detected. The dashed white lines highlight the spatial position along 
the slice (2500 km from the underlying umbra) where EM and T e time se- 
ries were created. These values are displayed in the lower panel, where the 
solid line represents the emission measure, and the dashed line traces the tem- 
perature. Both time series are displayed in their native units, as used in the 
top and middle panels. A clear anti-correlation between the EM and T e time 
series is readily apparent. 

conduction and optically th in radiation have been included 
dOwen & De Moortell l2Q09h . The ID model setup is con- 
structed to closely match the observed AIA 171 A emission 
at the lowest coronal part of the fan structure. In particular, 
the density profile (upper-left panel of Figure [9]) along the 
fan includes gravitational stratification, where the initial (z = 
0) value is determined through equalisation of the forward- 
modelled and observed 171 A intensities, as can be seen in 
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FIG. 9. — The initial (equilibrium) model density (upper-left) used in our ID LARE simulations, displayed as a function of distance along the coronal fan. The 
forward-modelled AIA 171 A emission (in DNpixel -1 ; solid line) is plotted in the upper-right panel, alongside the actual observed intensity profile (dashed line). 
Velocity amplitudes, induced by driving a slow magneto-acoustic wave into the numerical domain, are displayed in the lower-left panel, while the lower-right 
panel reveals the corresponding relative density perturbations, generated through subtraction of the forward-modelled values at t = 30 minutes from those at the 
beginning of the simulation. Damping, caused primarily by thermal conduction, is readily apparent in the lower two panels. All horizontal axes are in units of 
Mm, where 1 Mm = 1000 km. 



the upper-right panel of Figure [9] Here, the solid line corre- 
sponds to the spatial variation in 171 A emission, obtained 
through forward-modelling of the temperature and density 
profiles inferred from the observed fan structure (where the 
forward modelling is undertak en using a code developed by 
iDe Moortel & B radshaw 2008). A model temperature profile 
is constructed which utilises the minimum and maximum val- 
ues of the (isothermal) temperatures extracted directly from 
the AIA emission. The background magnetic field is taken to 
be constant. From the upper-right panel of Figure|9] it is clear 
that the theoretical model is in close agreement with the ob- 
served 171 A emission during the first 15 Mm along the fan. 
After this, the model emission declines more steeply than the 
observed 171 A emission. However, we are primarily con- 
cerned with the behaviour within the first 15 Mm, and there- 
fore can neglect larger distances where the emission becomes 
faint, and the observations become more affected by detector 
noise. 

A slow magneto-acoustic wave is driven into the domain 
at the lower boundary, through periodic perturbations of the 
field-aligned velocity component. A snapshot of the veloc- 
ity perturbations at t = 30 minutes, using a driving period of 
172 s to match that of the observations, is shown in the lower- 
left panel of Figure [9] A slow wave will propagate along the 
field at the characteristic tube speed, which is of the order of 
the local sound speed. Incorporating a temperature wl MK, 
and a background magnetic field strength salO G, the tube 
speed is of the order of 122 kmr 1 . Combined with a driv- 
ing period of 172 s, we expect a wavelength of the order of 
22 Mm, which can be verified in the lower panels of Figure[9] 



After an initial increase, due to the rapidly decreasing back- 
ground density, the velocity amplitudes are strongly damped 
through a combination of thermal conduction, optically thin 
radiation, and compressive viscosity. The corresponding den- 
sity perturbations (i.e., the density profile at t = has been 
subtracted from the density at t = 30 minutes) are shown 
in the lower-right panel of Figure [9] Such strong damping 
is predominantly caus ed by thermal conduction, as was pre- 
viously uncovered by [De Moortel & Hood! d2003l I2004D and 
lOwen & De Moortell (120091) . Only a small fraction of the 
overall damping can be attributed to optically thin radiation 
and compressive viscosity. 

To achieve a more meaningful comparison with the per- 
turbations observed by the AIA instrument, we utilised the 
model density and energy (temperature) distributions along 
the loop, to forward-model the emission for different AIA 
channels, at a multitude of time steps. We find that the mod- 
elled 171 A emission is significantly more intense over the 
first 15 Mm of the coronal fan structure, when compared to 
the 193A and 21lA intensities. Although these AIA chan- 
nels contain some emission lines at cooler temperatures, the 
193 A and 211 A bandpasses mainly r espond to hotter coro- 
nal temperatures ( i.e., above 1.5 MK; lO'Dwver et ail 120 lot 
iBrooks et alj|201 II) . and hence only begin to display elevated 
intensities when the density reduces and the local tempera- 
ture increases beyond 1 MK. This is consistent with the AIA 
observations, whereby the cooler 171 A emission dominates 
the first 15 Mm from the underlying sunspot. When forward- 
modelling intensity oscillations, it is important to remember 
that a variation in amplitude of the intensity perturbations is a 
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FIG. 10. — Time-averaged global wavelet power spectra, for the observed 
(top) and simulated (bottom) AIA time series, plotted as a function of distance 
along the coronal fan. The red, green, and blue lines highlight the oscillatory 
power contained within the 171 A, 193 A, and 211 A bandpasses, respec- 
tively, with each power spectrum normalised for clarity. Dashed green lines 
displayed in the observed spectra (top) relate to the power associated with 
the hot (~1.6 MK) and cool (<1 MK) spectral components of the 193 A 
channel, with the cool component closely following the 171 A trend, and the 
hot component tending more towards the 211 A spectrum. The damping of 
magneto-acoustic wave power, at hotter temperatures first, is an indication 
that thermal conduction is the dominant damping mechanism. 



combination of the change in amplitude of the model den- 
sity and temperature fluctuations, as well as the individual 
response functions of the AIA channels to different temper- 
atures. 

Finally, we utilise the global wavelet power, at each spa- 
tial position along the coronal fan, to compare the ampli- 
tude decay rates found in the different AIA bandpasses in 
a more quantitative way. To do this, a wavelet transform 
was performed at each position along the fan, with the 
temporally-averaged global wavelet power subsequently cal- 
culated. The red, green, and blue lines in Figure [10] repre- 
sent the time-averaged oscillatory power around 3 minutes, 
for the AIA 171 A, 193 A, and 211 A observed and forward- 
modelled bandpasses, respectively. Note that each of these 
have been normalised to their individual maxima for display 
purposes. Of particular note is the observed 193 A power 
spectrum. Here, the solid green line represents the oscilla- 
tory power registered through the entire AIA bandpass, while 
the dashed green lines relate to the power associated with 
the hot (^1.6 MK) and cool (<1 MK) spectral components 
of this channel. These comp onents have been separated us- 
ing the methods detailed in Kiddie et al] d2012l) . with the re- 
sulting cool emission peaking alongside the 171 A power at 
«4.5 Mm, and the hot emission displaying a peak between 
the 171 A and 211 A bandpasses (?«3 Mm). In reality, the 
global wavelet power of the 193Aand21lA emission is sub- 
stantially lower than that of the 171 A emission. However, 
the time-averaged power of the 193 A and 211 A modelled 
intensities are very similar, as is the case for the observa- 
tions. For these wavelengths, the simulated global wavelet 
power reaches a minimum around 10 Mm, while the 171 A 
bandpass reaches its minimum power slightly further along 



the fan structure, at about 15 Mm. These results agree qual- 
itatively with the observed time-averaged power, where the 
193 A and 211 A channels near-simultaneously reach a min- 
imum at about 9 Mm along the coronal fan, with the 171 A 
power diminishing completely by 15 Mm. This reiterates the 
importance of thermal conduction in the damping of coro- 
nal slow-mode waves. The increase in time-averaged global 
wavelet power, after wlO Mm in both the simulated and ob- 
served 193 A and 211 A time series, is due to the fact that, 
despite the decreasing amplitudes of the model density and 
temperature, the percentage intensity perturbations in these 
channels are actually increasing. This is a result of the in- 
creased sensitivity of the AIA 193 A and 211 A bandpasses 
to the > 1 MK temperatures which have been reached by this 
distance along the fan structure. The increase is not as appar- 
ent in the observed 193 A and 211 A time-averaged wavelet 
power spectra, which may be due to the decreased signal-to- 
noise levels found in these channels, especially when com- 
pared to the 171 A bandpass. 

4. OVERVIEW AND CONCLUDING REMARKS 

Here we present high-cadence observations of the solar 
atmosphere, obtained using the latest ground- and space- 
based facilities. Prominent oscillatory behaviour is de- 
tected throughout the optical and EUV image sequences, 
with remarkable similarities found between the detected wave 
modes. First, a number of UD structures in the photo- 
spheric umbra are found to display intensity oscillations with 
a «3 minute periodicity. These oscillations exhibit consider- 
able power, with regions encompassing the UDs displaying 
more than three orders-of-magnitude stronger power than the 
background umbra (Figure|3]l. Next, chromospheric intensity 
and velocity measurements were analysed for the presence 
of co-spatial and co-temporal oscillations. Such phenomena 
were detected, both on larger spatial scales, and with small 
central offsets with respect to the underlying photospheric 
oscillations. By considering the extension of magnetic flux 
tubes from the solar surface out into the upper solar atmo- 
sphere, a geometric expansion of only 76% in the radial di- 
rection and an inclination angle <40° allows the observed os- 
cillations to be interpreted as originating from within the same 
magnetic flux tube. 

Following examination of the phase lag between chromo- 
spheric velocity and intensity components, a V — I phase 
angle of —87 ± 8° was derived, allowing these waves to be 
described as a magneto-acoustic mode, with characteristics 
consistent with standing acoustic modes (middle panel of Fig- 
ure|5). The generation of a chromospheric standing wave may 
be the result of partial wave reflection at the transition re- 
gion boundary. An average blueshift velocity of wl.5 km s _1 
was found in the locations where high chromospheric os- 
cillatory power was present. This is a sub-sonic velocity, 
and coupled with a time lag between photospheric and chro- 
mospheric oscillatory phenomena, strengthens our interpreta- 
tion that the observed oscillations are upwardly-propagating 
magneto-acoustic waves, which originate in UD structures lo- 
cated in the photospheric umbra. 

A prominent fan structure is present in the simultaneous 
coronal EUV images, namely those from the AIA 131 A, 
171 A, 193 A, 211 A, and 335 A bandpasses. This fan is 
not readily apparent in either the transition-region dominated 
304 A emission, or in the higher temperature (~7.0 MK) 94 A 
bandpass. Using DEM techniques, we constrained the tem- 
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perature of the coronal fan to 0.5 - 1.2 MK, thus placing 
it outside the temperature range of both the 94 A and 304 A 
filtergrams (Figure [6). Time-distance techniques were em- 
ployed on the EUV imaging data where the fan structure 
was readily apparent, allowing the characteristics of propa- 
gating wave phenomena to be uncovered. Most coronal chan- 
nels, regardless of their absolute temperature sensitivity, re- 
vealed outwardly propagating wave fronts with an average 
periodicity and velocity of 172 ± 17 s and 45 ± 7 kms -1 , 
respectively (lower panel of Figure |7). The out-of -phase 
nature between the derived temperature and emission mea- 
sure signals indicates the presence of a compressive wave 
mode (Figure [8}. This, coupled with a sub-sonic wave speed 
(«45 kms -1 ), highlights the fact that these coronal phe- 
nomena are best described as upwardly-propagating magneto- 
acoustic slow mode waves. 

Employing numerical simulations, we were able to ac- 
curately simulate the behaviour of the coronal EUV emis- 
sion. Utilising input parameters derived directly from the 
AIA observations, forward-modelling techniques allowed us 
to evolve velocity, density, and emissivity values forward in 
time, creating a time series which could be directly compared 
to the AIA observations. Crucially, our simulations revealed 
that thermal conduction is the primary damping mechanism 
behind the dissipation of magneto-acoustic slow-mode waves 
in the corona. Other mechanisms, including optically thin ra- 
diation and compressive viscosity, play a secondary role in the 
damping of these oscillations. 

The fan structure observed in the AIA images, which dis- 
plays signatures of propagating magneto-acoustic waves, ap- 



Ahmad. I. A., & Withbroe, G. L., 1977, Sol. Phys., 53, 397 
Al, N., Bendlin, C, & Kneer, R, 1998, A&A, 336, 743 
Arber, T. D., Longbottom, A. W., Gerrard, C. L., & Milne, A. M., 2001, JCP, 
171, 151 

Aschwanden, M. J., & Boerner, P., 2011, ApJ, 732, 81 
Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D., 1999, 
ApJ, 520, 880 

Aschwanden, M. J., de Pontieu, B., Schrijver, C. J., & Title, A. M., 2002, 

Sol. Phys., 206, 99 
Aschwanden, M. J., Wiilser, J.-P, Nitta, N. V., & Lemen, J. R., 2008, ApJ, 

679, 827 

Aschwanden, M. J., Wuelser, J.-P, Nitta, N. V., Lemen, J. R., & Sandman, 

A., 2009, ApJ, 695, 12 
Aschwanden, M. J., Boerner, P., Schrijver, C. J., & Malanushenko, A., 2011, 

Sol. Phys., 384 
Ballai, I., Jess, D. B., & Douglas, M. 201 1, A&A, 534, A13 
Balthasar, H., Wiehr, E., & Kueveler, G., 1987, Sol. Phys., 112, 37 
Banerjee, D., O'Shea, E., Goossens, M., Doyle, J. G., & Poedts, S., 2002, 

A&A, 395, 263 
Beebe, H. A., & Johnson, H. R., 1969, Sol. Phys., 10, 79 
Bogdan, T. J., & Judge, P. G., 2006, Royal Society of London Philosophical 

Transactions Series A, 364, 313 
Bohlin, J. D., Purcell, J. D., Sheeley, N. R., Jr., & Tousey, R., 1975, BAAS, 

7, 356 

Brooks, D. H., Warren, H. P., & Young, P. R., 2011, ApJ, 730, 85 
Cauzzi, G., Reardon, K. P., Uitenbroek, H., et al., 2008, A&A, 480, 515 
Cavallini, F., 2006, Sol. Phys., 236, 415 

Centeno, R., Collados, M., & Trujillo Bueno, J., 2006, ApJ, 640, 1153 
Chou, D.-Y., Yang, M.-H., Zhao, H., Liang, Z.-C, & Sun, M.-T., 2009, ApJ, 

706, 909 
Cram, L. E., 1978, A&A, 70, 345 

De Moortel, I., 2006, Royal Society of London Philosophical Transactions 

Series A, 364, 461 
De Moortel, I., & Bradshaw, S.J., 2008, Sol. Phys., 252, 101 
De Moortel, I., & Hood, A. W., 2003, A&A, 408, 755 
De Moortel, I., & Hood, A. W., 2004, A&A, 415, 705 



pears to have anchor points in the south-west quadrant of 
the photospheric sunspot umbra (lower panels of Figure |2}. 
These locations are also consistent with the presence of 
large-amplitude wave phenomena detected in simultaneous 
photospheric and chromospheric image sequences. The co- 
temporal and co-spatial relationship between these upwardly- 
propagating magneto-acoustic wave modes, detected through- 
out the entire solar atmosphere, suggests such coronal phe- 
nomena may be driven by UD oscillations occurring inside 
the sunspot umbra. With this conclusion, it appears that pho- 
tospheric structures which are on the order of '.'5 (360 km) in 
diameter, can have a strong influence on coronal structures not 
only several thousand km above their position, but on struc- 
tures which have expanded into the local plasma to diameters 
often exceeding 10" (7000 km). 
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